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Abstract 

We discuss the possibility of implementing a universal quantum XOR 
gate by using two coupled quantum dots subject to external magnetic 
fields that are parallel and slightly different. We consider this system 
in two different field configurations. In the first case, parallel external 
fields with the intensity difference at each spin being proportional to the 
time-dependent interaction between the spins. A general exact solution 
describing this system is presented and analyzed to adjust field param- 
eters. Then we consider parallel fields with intensity difference at each 
spin being constant and the interaction between the spins switching on 
and off adiabatically. In both cases we adjust characteristics of the ex- 
ternal fields (their intensities and duration) in order to have the parallel 
pulse adequate for constructing the XOR gate. In order to provide a com- 
plete theoretical description of all the cases, we derive relations between 
the spin interaction, the inter-dot distance, and the external field. 



1 Introduction 

At present there exists a belief that single quantum systems, e.g., atomic traps, 
and electronic devices whose operation involves only a few number of electrons, 
could be used to implement the so-called quantum gates for quantum compu- 
tations. Quantum computations (quantum computers, or simply QC in what 
follows) can efficiently solve problems that are considered intractable by classical 
computers, e.g., the factoring of primes [1] and the simulation of others quan- 
tum systems [2]. As for classical computers, for the QC the accomplishment of 
an arbitrary algorithm can be performed using just a few specific manipulations 
called universal quantum gates. With these universal quantum gates, a process 
that acts in an arbitrary number of quantum bits (qubits) can be constructed 
using gates that act only in one and two qubits. 
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Although the first realizations of quantum computers were performed using 
trapped ions [3], NMR [4], and optical cavities [5], it is believed that, as for 
classical computers, for QC the most promising object for a possible large-scale 
implementation, i.e., a implementation involving a large number of qubits, will 
be the solid-state devices. Among these devices, one can highlight the system 
of two coupled semiconductor quantum dots (QD) [6J. In this system the qubits 
are the one-half spin states of an excess electron in each dot, and the universal 
quantum gates can be obtained by performing arbitrary rotation of the individ- 
ual spin (one-qubit gate) and one operation capable of entangling the spins of 
the two electrons (two-qubits gate). Any such operation with this entangling 
feature can be used [7j- The rotations and the entanglement operations can 
be performed by applying external electromagnetic pulses to the dots. These 
pulses can be applied in a serial manner, which perform sequentially the one and 
the two-qubit gates using a sequence of pulses, or in a parallel manner, where 
the two and one-qubit gates are performed at once by the application of just 
one adequate pulse (called a parallel pulse). Due to the decohering effects it is 
important to perform the operations in a minimal time, which, as discussed in 
PD, can be achieved more adequately with a parallel pulse. 

The adequate parallel pulse to implement a quantum gate can be obtained if 
we know the Hamiltonian that describes the coupled QD, but to construct this 
Hamiltonian, one has to know the exact relation between the external param- 
eters of the system (such as the external fields and the physical characteristics 
of the dot) and the interaction between the electrons. An advance in the de- 
scription of this relation was obtained in [9], where the authors have used the 
Heitler-London approximation to obtain an explicit expression that relates the 
interaction between the dots, the inter-dots distance and the external electro- 
magnetic field. However, in this article, the two dots were subject to the same 
magnetic field, and, as we will see, a little difference between the fields in the 
dots possess a great influence in the evolution of the system and, consequently, 
in the implementation of the gates. In addition, after the Hamiltonian is known, 
we have to construct exact solutions of the evolution equation of this four-level 
system. In our work [10] we studied possible exact solutions of a system of two 
coupled spin one-half particles subject to external time-dependent magnetic 
fields. 

In the present work we discuss one possibility of implementing an univer- 
sal quantum gate, namely the exclusive or (XOR) gate, by using two coupled 
QD subjected to external magnetic fields that are parallel and slightly different. 
First, in Sections[2]and[3j we make a brief discussion of the model with the paral- 
lel external magnetic fields of general form, and describe possible corresponding 
exact solutions from our previous work [TO]. In order to have a complete theo- 
retical description of the problem, one needs to find relations between the spin 
interaction and external fields. This problem is solved in Section |H To this end 
we generalize results of [9], derived for equal parallel external fields, to the case 
of different fields, obtaining relations between the spin interaction, the inter-dot 
distance, and the external field. In section [5] we consider two parallel fields with 
specific form. In the first case, parallel external fields with intensity difference 
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at each spin being proportional to the interaction between the spins, the latter 
depending on time. General exact solution describing this system is presented 
and analyzed to adjust field parameters. Then we consider parallel fields with 
intensity difference at each spin being constant and interaction between the 
spins switching on and off adiabatically. In both cases we adjust characteristics 
of the external fields (their intensities and duration) in order to have the parallel 
pulse adequate for constructing the XOR gate. 

2 Four-level system 

In non relativistic quantum mechanics the dynamics of a fixed spin one-half 
particle, subject to a time-dependent external field K(f), can be described by 
the Hamiltonian h = (<r • K) [11], where cr = (ai, <T2, C3) are the Pauli matrices 
and we set h — 1. The quantum Hamiltonian of two interacting one-half spins 
subjected to external fields G and F, respectively, is chosen as, see e.g. [12], 

ff(G,F,J) = (p.G) + (£.F) + |(£.p) , (1) 

where the function J — J (t) defines the spin interaction, G = {G\ (t) , G2 [t) , G3 (t)) 
and F = (Fi (t) , F2 (t) ,Fs (t)) are the time-dependent external fields at each 
particle, and the 4x4 Dirac-matrices p and £ are in the standard representation 

3 

£ = /<X>er,p = er<X)/,(£-p)=er®£r = ^cr i <X>cr 4 , (2) 

where / is the 2x2 identity. 

The interaction of the spins in {1} is known as the Heisenberg interaction 
and describes, for example, the interactions in the well known Hubbard model 
[13] , In particular, under certain conditions, this interaction can be used to 
describe the coupling between two QD [6]. In this case, it is possible to control 
not only the external fields, but also the spin interaction. When these QD are 
used to implement a quantum gate, the operation of the gate is performed by 
varying the external fields and the spin interaction during a certain time r. To 
this end, in order to justify the Heisenberg interaction, the following conditions 
(see [6]) must hold: 

1. The time r cannot be too small to avoid transitions to higher energy 
levels, so that the time scale should be bigger than H/AE, where AE is 
the difference between the first two energy levels; 

2. The decoherence time of the physical system should be much bigger than 

T. 

The use of these QD to implement quantum algorithms requires controlling 
the individual spins, which can be done by the fields F and G, and any inter- 
action between the spins capable of creating a entangled state starting from a 



3 



original product state [7]. A system described by the Hamiltonian (TJ satis- 
fies these conditions. Namely, when the fields F and G are zero, the evolution 
operation R t (G, P, J) of the Hamiltonian ((J) can be written as [10] 

R t (0, 0, J) = exp [i$ (t) /2] [I cos $ (t) - ^ sin $ (t)] , 
A=l[I+(S.p)] , $(*) = [ J(r)dr, (3) 

z -/to 

where I is the 4x4 unity matrix. From the above expression we see that, when 
$ = 7r/4, the evolution operator acts as the universal gate known as square root 

1/2 

of swap (U S w )• The interaction function J can experimentally be controlled in 
many different ways, e.g., by applying electrical or magnetic fields to the dots 
|14| . So we can construct any quantum gate by a sequence of pulses by turning 
on and off the external fields and the interaction. For example, the XOR gate 
can be constructed as [9] 

Uxor = exp (i Jp 3 ) exp (-«|s 3 ) U^ 2 exp (i^Ps) U^ 2 . (4) 

On the other hand, the construction of gates by a sequence of pulses is not 
appropriate, because the duration of all the sequence can be too long, violating 
condition (2) above, or the pulses need to vary too fast, violating the condition 
(1). So it is important to be able to implement the gates at once, applying just 
one adequate field. This single pulse is called parallel pulse [§]. For example, in 
the case of the XOR gate (@]), to use a parallel pulse one needs to find a field 
whose evolution operation, at a given instant r of time, has the form 

i? T (G,F,J) = [/xoR = exp[-zJ(£ 3/93 + X 3 + /93)] . (5) 

In order to find this parallel pulse, in a general case, we need to know exact 
solutions of the Schrodinger equation with the Hamiltonian (TJ for different 
kinds of external fields and interactions. A large number of these exact solutions 
are present in our previous work [10], and here we will use some of these results 
to describe the implementation of a £/xor gate. The procedure developed here 
can easily be extended to others gates using others fields present in |10| . 



3 Parallel external fields 

Let each spin in our system be subject to different time-dependent external 
magnetic fields along the z direction, 

G = (0, 0, nagiBi) , F = (0, 0, MflfiaBa) , Bi,a = B h2 (t) • (6) 

where \ib is the Bohr magneton and the g-factor of the dot i. Therefore the 
Hamiltonian ([1]) assumes the form 

H = i [(S 3 + pa) B+ - (S 3 - p 3 ) B- - J] + AJ , B ± = fx B {g x B x ± g 2 B 2 ) , 

(7) 
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with the constant 4x4 orthogonal matrix A given in (J3j) . For the first and 
fourth components V\ and of the four-spinor '5, solution of the Schrodinger 
equation (i^ = H^) with the above Hamiltonian, we get 



Vi = C\ exp 



B+ dr 



, V4 = C4 exp 



- B+ dr 



(8) 

where are complex constants. Besides, for the components ^2,3 of we 
obtain the expression 



itp' = 



(cr-K) — ^ 



''2 

t>3 



K(i) = (J(t),0,B_ (t)) 
Making the transformation 



t// (t) = exp 



(9) 
(10) 

(11) 



the two-component spinor ip obeys the equation for a single particle with spin 
one-half in an effective external field K (t), i.e., iip = (cr-K) ip. Consequently, in 
this case, the four-level system problem reduces to finding solutions of a two- 
level system. Writing the solution of this two-level-system as if) (t) — u t tp{0), 
using the 2x2 evolution operator u t (J, B-), we can write the evolution operator 
of the four-level-system governed by the Hamiltonian ([7]) as 



(12) 



Rt (G, F, J) = exp \-- [(S 3 + pa) T (t) + £ 3 p 3 $ (t)] ) M (t) 

r(t)= / B+(t) dr, $(*) = / J(i) dr, 
Jo Jo 



with the 4x4 matrix M (t) given by components of the single particle evolution 
operator u t as 

10 
M = I u t 
1 



Writing the effective field K ([TO]) as 

q+ [q x q] 



K 



1 



(13) 



for an appropriate real three- vector q(i), we can obtain the evolution operator 
Ut in the general form [TT] 



ut = 



1 + qq - »gp 

v/(l + q 2 )(l + qj5) 



, p = q - q + [qo x qj 



(14) 
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where qo = q(0). Thus we can construct the XOR operator |(5]) by using any 
effective field K given by a vector q that, after a certain time T, points again 
in the initial direction (q (T) oc q ). When this is the case, ut — I and the 
evolution operator R t fT2|) will assume the form Uxor © when 

r (T) = $ (T) = | mod (2tt) . 

This periodicity in the direction of the vector q can be considered a general 
property for which the parallel field acts as a XOR gate. An arbitrary choice 
of this vector, respecting only this periodic condition, can be used to obtain a 
variety of parallel fields that will certainly possess the adequate characteristics. 

Although the interaction function J depends on the applied external fields, 
as we will see in the next section, when £?_ << B + , we can make J — J (B+) 
and consider the vector K in (fT0| as composed of two independent functions 
J (t) and B- (t). Besides, the interaction function can be controlled by electric 
fields [14], whose interference in the spin states via the spin-orbit coupling, in 
many practical application, can be neglected. Then, if necessary, we can consider 
the functions J, £?_ and B + independently. In addition, although the obtained 
expressions depend on the sum and on the difference between the fields at the 
spins, only the average value of this sum is relevant, so that the explicit form of 
its time variation can be arbitrary. 

Some comments on the experimental realization are in order. With the 
technology available nowadays, systems of two coupled QD, each holding ex- 
actly one electron, are routinely constructed in experimental solid state physics 
|15| . The confinement of the electron and the control of the interaction between 
them can be performed by electrical gates directly assembled in the semiconduc- 
tor material using etching techniques. The coupled electron spin is a result of 
the combination of the Coulomb interaction and the Pauli Exclusion Principle. 
The spin state of this system can be initialized using a variety of techniques 
[Hllll]. In all experiments with such systems, the difference between the en- 
ergy of singlet and triplet spin states are controlled by applying either static 
or time-dependent external magnetic fields of intensity up to the order of a few 
teslas. See, for example, the experimental setup [14]. Therefore, the construc- 
tion of systems of two-coupled QD subject to external parallel magnetic fields 
is commonplace. However, the main problem in the experimental application 
of the results presented in this work is to generate, and control, a difference 
in the magnetic field felt by each spin, that is, the manipulation of B- ([7]). 
The quantity £?_ is not necessarily associated with a gradient in the applied 
magnetic field. It is possible to couple two different QD [17] in such a way that 
5i 7^ 92- In addition, some techniques permit the manipulation of the g-factor 
by changing the size of the dots or by the application of external electromag- 
netic fields [IHIIIS]. The hyperfine coupling between the electron spin and the 
nuclear spins of the semiconductor material can also be explored to obtain a 
field gradient by producing a differential Overhauser field [20]. As a result, we 
can have a nonzero i3_ even when the two dots are subject to the same external 
magnetic field. Besides, it is possible to fine-tune the quantity B- with the 
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application of a localized magnetic field [21]. In ]22J an experimental setup of 
this kind is used to control the spin state of a single electron in a coupled QD 
system. The fine-tuning of the quantity B_ has been an intense object of study 
in the current solid state experiments. 

4 Relation between the interaction function and 
the external fields 



In order to construct the evolution operator (|T2|) for some specific parallel field, 
one needs to know the explicit dependence between the external fields and the 
interaction function. In the article ]9], the authors study two QD, coupled with 
the Heisenberg interaction (Q}, and use the Heitler-London approximation to 
obtain an expression for the interaction function J (V, a, B, E) as a function of 
the gate voltage between the coupled QD (V), the inter-dot distance (a), an 
external magnetic field (B) and an external electric field (E). The knowledge 
of the interaction as a function of various parameters allows to compensate the 
changes in J with the variation of one given parameter by controlling another 
parameter. So one can, for example, vary the magnetic field B and maintain J 
constant by changing the inter-dot potential V. This procedure allows treating 
some parameters of the Hamiltonian of the system as independent. However, 
the analysis developed in [9] concerns only the problem of two dots subject 
to the same magnetic field. Since the work developed here depends on the 
difference between the fields at each dot, we will repeat the process given in 
the cited article, but with different magnetic fields at each dot. For different 
parallel magnetic fields (i?i,B 2 ) applied at the spins, the non-homogeneity of 
the field make it very complicated to calculate J (Bi, B2), since, in this case, the 
Zeeman terms (p-G and X-F) will not be independent of the space coordinates 
anymore. Despite this fact, since only a small field difference is important, the 
fields can be considered homogeneous inside each dot and the Heitler-London 
approximation can still be used to describe the interaction function. 

The model consists of two identical bi-dimensional dots with an inter-dot 
separation of 2a, each dot with one electron of charge e and subject, respectively, 
to the magnetic field B\ and B2 in the z direction. The field Bi is considered 
homogeneous inside the dot i. In the Heitler-London approximation, we start 
with a combination of the orbital ground states of bi-dimensional single dots 
and combine this state in symmetric (| X I > +)) and anti-symmetric (|^-)) states 
for the double dot problem as 

iw = '- +)±l+ -> 



where <p± (r) = (r |±) denotes the one-particle orbital centered at r = (±a, 0) 
and S = (+ |— ) the overlap between the two orbitals. The interaction function, 
or the exchange energy, is the difference between the energies, J = (\?_ | H or b 
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(vE'+l H or b where the double-dot orbital Hamiltonian is given by 
H orb = h° (n) + h° (r 2 ) + C , C 



« ri - r 2 



ft (r<) 



1 

2m 



-A(r 4 ) + V(r ( ) , i = l,2 , 



A ( r i) = y (~S/t) ^i, 0) , 1/ (x, y) 



1 



4a 2 ^ 



(15) 



In this expression c is the speed of the light, k is the dielectric constant of the 
medium, m is the effective mass of the electron (e.g., m = 0,067m e in GaAs), 
r, and pi are the position and the momentum of the i-th electron, Bi is the 
magnetic field at dot i and the harmonic potential well V (r) of frequency loq is 
motivated by experimental results [23j • The matrix elements needed to calculate 
J can be obtained by adding and subtracting in 1|15|1 the harmonic potentials 
V± centered at Xi = (— l) 1 a for the i-th electron, 

H orb = h Q _ (n) + ft!}. (r 2 ) + W + C , 
4^) = ^ [pi-^A(r i )] 2 + V ± (r i ) , 



^± fa) 



(si±o) 2 + »?l , W = J2V(r i )-[V-(r 1 ) + V + (r 2 )\ , (16) 



i=l 



and using the ground states functions ip± centered at r = (±a,0), 



¥>± (r) 



exp ± 



mu>± 



y exp 



2/) 



w± = A /^5 , w L _ (+) = eB 1(2 )/2mc , 



which are eigenfunctions of /irj. with energy w± , 

h°_ (+) (ri( 2) ) (ri (2) ) = cj_ (+)¥ >_ (+) (r 1(2) ) 

With this, the expression for the interaction function becomes 



J = 



2S 2 



L 



htu (b\ - b 2 _) (&_ - ^ 



(1-5 4 ) 
L = ( (12\C + W\12) - 



4 b+b- 

Re (121 C + W 121) 



S' 2 



using the expression ifTTj) for the functions ip± (r^) = (r$ |±) we get 

S 2 1- A 2 

(1 - S 4 ) ~ (2 sinh (2Af ) + A exp (-2M) (2 - A 3 )) ' 



M 



2d 2 
6+ + bi 



b-b+ + 



4cj 2 



b- - hj 
6- +6h 



n 
«o 



(17) 



(18) 
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where ao = ^h/mujo is the effective Bohr radius of the dot. For the matrix 
elements of W in ifTSl) we have 



(12| W\12) 
2 



(121 W 121) 



S' 2 



2d 2 (b- + ^ 



1 + A 2 
(1-A 2 f 



1 



A 2 - 1 \ d / a 4 .2 \ 

i r- ) — 7t (A — 6A 2 - 3) 

b- + b+J 2 v ; 



and for the matrix elements of the electric interaction C between the electrons 
(12| C |12> ««M<7|21) 

,2 



s 2 



— J^b V / (1-A 2 )exp [-d 2 (1 - A 2 ) 6] I [d 2 (1 - A 2 ) b\ - exp 











Io 




2 


2 



^ = 6(l + A 2 )-i + y[(l-A 2 )6] 2 -2(l + A 2 ) + ^, b= h -±^±, 

where Io is the zeroth-order Bessel function. The quantity A is related to the 
difference and to the sum of the fields (B'± — B\ ± B2, without the g-factor) at 
the dots by 



2h^- 

MB 



B 



1 2 



B 



1 2 



2h^- 

fJ-B 



-1 2 



B 



1 2 



B 



1 2 



- (2B'B'_ 



This quantity can be used to evaluate the error in considering J independent 
of B'_. In the denominator of the above expression, usually only the first term 
is relevant, e.g., for a typical GaAs dot we have Swo/mb = 50 T. For the same 
field at both dots we obtain the result given in [9]. 



5 The exclusive OR gate 

Here we establish the parameters of two kinds of parallel external fields that can 
act as a XOR gate, that is, fields in the form J6]) whose evolution operator (fT2|) . 
at a given instant of time, has the form J5]). The same analysis presented here 
can also be carried out for the 26 different families of external fields presented in 
[TT] , which gives a wide range of fields to choose from that are more appropriate 
to adjust to experimental setups. All the development depends on the ability 
to vary independently the difference i?_ and the interaction J, but, as we saw 
in the previous sections, there are various circumstances when this requirement 
can be satisfied. 

5.1 B- proportional to the interaction 

Suppose that we obtain, theoretically or experimentally, an expression for the 
interaction function J, and we choose B~ such that it is proportional to J. As 
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a consequence we can write 



J (t) = q (t) sin A , B-(t)=q (t) cos A , 

where q (t) is an arbitrary function of time, and A is a real constant. Then, 
the general solution of equation J9]) can be written as i\> (t) = utip (0) with the 
evolution operator ut given by 



u t = cos u> — i (<7i sin A + 03 cos A) sin ui , li (t) = / q (t) dr 

Jo 



(19) 



Selecting the instants T when cos a; = 1, 



u(T) = I q (t) dr = 2rm , n e N* , (20) 



the evolution operator (fT2| of this problem assumes the form 
Rt = exp (i3nir) exp 



2 

2 ( r I s 3 + P3] + 2n7rS 3 p 3 sin A) 



So to obtain the XOR gate J5]), up to a phase, at every instant T i[20|) . we have 
to set 

T (T) = 2nvr sin A = ^ mod (2tt) , 
or, in a more explicit form, 

•a 4m + 1 

sin A = , m s -N , m < n . 

in 

When the interaction J and the difference £?_ are constants, so q is a con- 
stant, we obtain 

r(T) = jrmod(2^) , T= . (21) 

I J 2 + B 2 _ 



For the special case of a constant sum of fields B + , the above expression can 
be compared with the XOR gate obtained in [8]. To have a idea about the 
physical values involved in the expression lj2lj) . for a typical GaAs QD [9], for 
an interaction about 50 /ieV and a £>_ of 10 mT we obtain a XOR gate in a 
time of 10 ps. 

5.2 Adiabatic pulse 

The implementation of the quantum gate is performed by the variation of the ex- 
ternal fields. However, as described in condition (1) of section[2l these variations 
can not be so fast so as to prevent the excitation of higher energy- levels. This 
problem can be avoided by using an adiabatic variation of the fields, which can 
be obtained by a time-dependence in the form of sech (ut) [8], with oj « AE/h 
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(see condition (1) in section In this section we analyze the case when the 
field's difference is constant and the interaction varies adiabatically. The case 
when both the interaction J and the difference B- vary adiabatically is a spe- 
cial case of the preceding section. The case when the interaction function is 
maintained constant and the quantity -B_ vary adiabatically can be study by 
knowing that [TT] : if ip is a solution of the equation ([9]) with the external field 
(K 1 ,0,K 3 ),then^> = (2)- 1/2 (a 1+ a 3 )^ is a solution of |(9]) with the external 
field (-K3, 0, K\). In addition, for £?_ not proportional to J a lot of solutions of 
this sech form can be found in [IT] . 

Suppose that, by vary B + or by controlling the electric potential between 
the dots, we obtain a variation in the form 

J(t) = a/coshut , B_ = c , (22) 

where a, c are constants. The evolution operator u t of equation ([9]) with the 
above fields is given by 



|G°| 2 + |G?| 2 V G 2 ( Z ) G|(z) J { -G° G\ 

where the * indicates complex conjugation and 

Gi (z) = (2c - iw) aT" (1 - z)" F (A, -A; 7; z) , 

G 2 (z) = 2az-» +1 ' 2 (1 - z)" +1/2 F (1 + A, 1 - A, 7 + 1; z) , 

z (*) = ~ (1 - tanhwt) , G° = Gi (1/2) , 

7=1/2-21/, A=H. (24 ) 

In these expressions, F (a, /3, 7, z) is the Gauss hypergeometric function. 

If we use the sech variation to create a pulse of the interaction function, this 
pulse will be turned off when t » 1/lu. In this time limit we have lim wt ^oo z = 
e -2wt an( j p ^ ^ ^ q-j _ ^ jf m this limit, we choose 

G° 2 = aF (l + A,l - A,7 + l;0 =0 , (25) 

the evolution operator (|23| assumes the form 

lit (u>i — » 00) = exp {—ia^tc) . (26) 

So, following the procedure of the preceding sections, we see that the evolution 
operator f(T2|) will behave as a XOR gate J5]) in the instant T, when T (T) = 
7r/2mod(27r) and, for c 7^ 0, 

lot: (1 + 4m) „ ?w 

T= — , n, m e N , (27) 



4 arctan [exp (wT)] — 7r ' c 
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or, in the limit ujT >> 1, a ~ u> (1 + 4m). 

The condition f25|) can be obtained in various ways, for example, in the 
special case when c = (and consequently lim^t^oo u = I in (|26|) and we can 
abandon the condition T = rar/c) we can use the relation [24] . 



So, we obtain the desired condition whenever \a\ — 2muj, m 6 N*. But, once 
lf27j) gives a as an odd multiple of u, it seems that the gate can not be obtained 
with equal fields at both dots. 

6 Some final remarks 

We have demonstrated that the general solution of the Schrodinger equation 
for the system of two interacting spins placed in parallel magnetic fields can be 
reduced to the general solution of the corresponding equation for only one spin 
in an effective external field. This allows us to use the known exact solutions 
of the latter problem to obtain exact solutions for the system of two spins. In 
turn, this was used by us to describe a construction of a universal quantum gate. 
We have found restrictions on external magnetic fields (in terms of the effective 
field for the one-spin system) such that the four level system can act as a XOR 
gate. The explicit form of the evolution operator and the known dependence of 
the spin interaction with the magnetic field show that a small difference in the 
fields applied at each dot (or a small difference in g-factors of these dots) can 
be explored to control the quantum gate without changing the spin interaction. 

One of the problems in the practical application of the theoretical results in 
constructing universal quantum gates is the knowledge of the spin interaction as 
a function of external fields. This function can be obtained, in a way unrelated 
to the Heisenberg interaction model, using a theory that takes into account the 
physical characteristics of the system, as we do here using the Heitler-London 
approximation. If from the expression of J, obtained from this more complex 
model, follows that J does not depend on B- for £>_ <C B+, we are free to 
choose £?_ as a function of J (B + ). By using this expression of J(B+), and 
choosing £?_ (<) to be proportional to J (B + ), we can measure the probability 
transition between the states |t|) and |J.|) (the swap operation). According to 
our results, this probability is equal to sin 2 A sin 2 w. This result can be used 
to test the consistence between the J(B+), obtained by any model, and the 
Heisenberg interaction model. 

Acknowledgement 1 M.C.B. thanks FAPESP; D.M.G. thanks FAPESP and 
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